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Abstract 

We introduce a simple representation for irreducible spherical tensor opera- 
tors of the rotation group of arbitrary integer or half integer rank and use these 
tensor operators to construct matrix product states corresponding to all the vari- 
ety of valence-bond states proposed in the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
construction. These include the fully dimerized states of arbitrary spins, with 
uniform or alternating patterns of spins, which are ground states of Hamiltonians 
with nearest and next-nearest neighbor interactions, and the partially dimerized 
or AKLT/VBS (Valence Bond Solid) states, which are constructed from them by 
projection. The latter states are translation-invariant ground states of Hamilto- 
nians with nearest- neighbor interactions. 
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1 Introduction 



The problem of introducing exactly solvable models in quantum spin chains has a 
long history in statistical mechanics and mathematical physics, which encompasses 
a variety of models and techniques [Tj '2] . One can mention the XY [3] , the Heisen- 
berg XXX and XX Z [4j, the AKLT [5j, and the Majumdar- Ghosh models [B], and 
the free fermion, the Bethe ansatz, and the matrix product techniques to name only 
a few of the most important models and techniques which have been developed so 
far. Finding any new exactly solvable model, is an important step, since it acts as a 
reference model for developing approximate perturbative solutions for more realistic 
models. It will also help us to test many of the new ideas about collective behavior 
of quantum systems, i.e. entanglement properties [3 [8], or the relation of criticality 
and universality of entanglement [9l [TO] . 

In [^, Affleck, Kennedy, Lieb and Tasaki (AKLT), suggested a new construction 
for a variety of spin states, known as valence bond states. The basic element of this 
construction is a spin-1/2 singlet state, a dimer, \s) = ~ which is 

called a valence bond in [5j. A dimerized state is just a juxtaposition of such dimers 
on adjacent sites, figure ([T]-a). Such a state is clearly seen to be a ground state of 
a Hamiltonian with three-sites interactions (nearest and next-nearest neighbors), the 
local Hamiltonian of which is the projector to spin 3/2 states, h = The reason 

is that due to the presence of a dimer, the sum of spins of three adjacent sites adds 
up only to spin 1/2. The parent Hamiltonian of this fully dimerized state, is known 
as Majumdar-Ghosh Hamiltonian and has the form 



This Hamiltonian has a two-fold ground state degeneracy, the other ground state be- 
ing simply a one-site translation of dimers to the left or right. 

One can also consider fully dimerized states [5] with alternating patterns of spins, 
where there are alternating number of valence bonds or dimers. An example of this 
is shown in figure ([2]-a), where the local three-sites Hamiltonian, should be taken as 
projector to spin 2, h = P2- Moreover one can use projection, to construct from these 
fully dimerized states, partially dimerized or AKLT/VBS states which are ground 
states of Hamiltonians with nearest-neighbor interactions. For example in figure ([T]- 
b), if one projects each pair of spin-1/2 particles in a bulb of the original chain to the 
symmetrized triplet, a non-dimerized spin-1 state is obtained on a new chain, whose 
parent Hamiltonian which annihilates this state is the sum of spin-2 projectors P2 on 
consecutive sites. The reason for this annihilation is that the sum of four initial spins 
on the original chain (known also as the virtual chain) add up to at most spin 1, due 
to the presence of the valence bond which is a singlet. In this way a spin-1 quantum 
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Figure 1: Color Online. A fully dimerized valence bond state, and the AKLT con- 
struction of a partially dimerized spin-1 state. The states in each bulb are projected 
to the symmetrized triplet. The parent Hamiltonian of the upper chain (a) has an 
interaction range of 3 lattice sites, while that of the lower chain (b) has a range of 2 
sites. 

chain is obtained which is the exact ground state of the following Hamiltonian: 

H = Si ■ 5i+i + -{Si ■ (2) 
j 

Projection can also be used for other types of dimerized state as shown in figure 
([2]-b) to construct states with arbitrary integer or half integer spins. For example in 
figure dU-b), looking at the number of valence bonds which are singlets and are not 
counted in the addition of spins in the virtual sites, one finds that the local Hamilto- 
nians can be chosen as hi = X3P3 and /12 = A3P3 + A2-P2, where Pj is the projector 
on spin-j states and Aj's are positive coefficients. To assure translation invariance for 
the parent Hamiltonian one then takes H = where h = P3, is the operator 

common to both hi and /12. Needless to say, this construction can be generalized 
by taking different alternating number of dimers in the virtual chain. This is also 
the basic idea behind the exactly solvable spin-3/2 spin systems on the honeycomb 
lattice or more generally the basic idea behind PEPS, or Projected Entangled 
Pair States [12], which has only recently been discussed in the literature. 

In the course of time, the basic idea of AKLT, which in turn was inspired by 
the work of Majumdar and Ghosh led to the development of finitely correlated 
or matrix product representation of states [13 ^ 114 ^ [T5]. a representation which when 
existing, greatly facilitates the calculation of many properties of the ground states of 
quantum systems [gil5l[IS[I71lig[IS[Ml[2ll[22l[l3l[Ml[23|26]. 
The Matrix Product (MP) representation was also found to be closely related with 
the success of density matrix renormalization group |27l [28l [29] . 

When considering spin chains, the basic continuous symmetry is the rotation sym- 
metry captured by the su{2) group, and there has been many different and equivalent 
implementations of this symmetry in matrix product states |2 H 122 ^ 125 ] 1261129] . While 
a lot of progress has been made in defining matrix product states, having specific 
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Figure 2: Color Online. An example of a dimerized state with alternating spin pattern 
and the AKLT construction of spin 3/2 chain. The interaction ranges are similar to 
that of figure ([T]) . 

symmetries, to our knowledge the original AKLT variety of states, have not been cast 
into a simple and uniform matrix product form for both integer and non-integer spins. 
For the integer case however, such a formulation has been reported in [30]. There is 
no doubt that such a representation, will be of utmost importance for further study 
of AKLT models, and even for similar models on more general geometries like the 
Bethe Lattice [3T] . 

To ensure invariance of the Matrix Product State (MPS) under rotation, it is 
sufficient that the elementary matrices used in the definition of the MPS constitute 
a representation of spherical tensor operators of a specific rank. The rank of the 
tensor depends on the spin of the actual lattice and the dimension of the represen- 
tation determines the dimension of the auxiliary matrices. Finding a simple and 
minimal-dimensional representation for such tensors, constitute the basic problem in 
constructing rotationally invariant MPS, both for spin chain, spin ladders, or two 
dimensional lattices. 

What we will do in this paper is to provide a uniform and simple matrix product 
representation for all the AKLT, or valence bond states and even more general states. 
The starting point of our analysis is a simple and compact representation of spheri- 
cal tensor operators of any rank, integer or half integer. These tensors enable us to 
define MP representations for Majumdar-Ghosh states, (which are the ancestors of 
AKLT states) and their generalization to arbitrary spins, and then we will use them 
to construct MP representation for partially dimerized states. We then use projection 
method to find MP representations for arbitrary spin chains, with nearest neighbor 
interaction. The parallel with the AKLT construction is simple: the basic idea is to 
replace a collection of 2s spin-1/2 dimers or valence bonds with a single spin-s valence 
bond and represent the states constructed from these spin-s valence bonds as MPS. 

Besides having the benefit of calculability, when we have an MPS representation, 
the very method of MPS allows us to find a larger family of Hamiltonians than the 
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AKLT method. This larger family, with its larger number of couplings will enable us 
to better adjust or approximate an exactly solvable Hamiltonian with realistic situa- 
tions. We will see an example of this in this paper. 

The structure of this paper is as follows: In section[2]we review the matrix product 
formalism |13[ [T^ [T5] in a language which we find convenient [21] for further devel- 
opments. In particular we emphasize the symmetry properties of the ground state 
and the Hamiltonian. In section [3] we will introduce a compact formula for spherical 
tensors of rank s (integer or half integer) and use it to construct dimerized states of 
arbitrary integer or half-integer spins in section These are the generalization of 
Majumdar-Ghosh states, or fully dimerized states, to arbitrary spins. We then go on 
in section ([5]) to define MP representations for other types of dimerized states. In sec- 
tion ([6]) we find MP representations for AKLT/VBS states. The core of this section is 
the definition of new kinds of tensors, which play the role of auxiliary matrices for the 
MP representations of these states. Section ([7|) is devoted to some specific examples, 
where more detailed properties of some of the states and their parent Hamiltonians 
are derived. We conclude the paper with a discussion. 



2 Matrix Product States 

Let us first make a quick review of the matrix product states in a language which 
we find convenient [211 122| . For more detailed reviews of the subject, the reader can 
consult a more comprehensive review article like [23] or any of the many works where 
specific examples have been studied Ii[2ll[22l[ia[IZl[ia[2Ql[l9l[2l[iai26j. 

Consider a ring of sites, where each site describes a d— level state. The Hilbert 
space of each site is spanned by the basis vectors i = 0, • • • , d — 1. A state 

l^>= E ^hi2-iJki2---iN) (3) 

is called a matrix product state if there exist D dimensional matrices Ai G C^^^, i = 
0, • • • , d — 1, such that 

where Z is a normalization constant. This constant is given by Z = tr{E^), where 
E = Ef=o 4* ® 4- Note that we are here considering homogeneous matrix product 
states where the matrices depend on the value of the spin at each site and not on the 
site itself. More general MPS's can be defined where the matrices depend also on the 
position of the sites [23j . 

The collection of matrices {Ai} and {fiUAiU~^}, where // is an arbitrary complex 
number, both lead to the same matrix product state, the freedom in scaling with //, is 
due to its cancelation with Z in the denominator of ([3]). This freedom will be useful 
when we discuss symmetries. There has been discussions on the symmetry of matrix 
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product states in the literature [131 El US Ell ESI ESj , here we use the language or 
notation used in [21]. 



2.1 Symmetries of the ground state 

Consider a local continuous symmetry operator R acting on a site as R\i) = Rji\j) 
where summation convention is being used. R is a d dimensional unitary representa- 
tion of the symmetry. A global symmetry operator TZ := i?®^ will then change this 
state to another matrix product state 

Kh-iN ■= t^iKA^ ■ ■ ■ ^'^M)' (5) 

where 

A'i := RijAj. (6) 

The state l^*) is invariant (i.e. a singlet) under this symmetry if there exist an operator 
U{R) such that 

RijAj = U'^{R)A,U{R). (7) 

Repeating this transformation and using the group multiplication of the transforma- 
tions i?, puts the constraint 

Ur>Ur = Ur.r. (8) 

Thus U{R) is a L> dimensional representation of the symmetry R. In case that R is 
a continuous symmetry with generators T^, equation ([7]), leads to 

{Ta)^JAj = [Ta,Ail (9) 

where Ta and 7^ are the d— and D— dimensional representations of the Lie algebra of 
the symmetry. 



2.2 Symmetries of the Hamiltonian: 

Given a matrix product state, the reduced density matrix of k— sites is given by 

triiA* ■■■A* ^ Ai, ■ ■ ■ AiJE^-'') 

The null-space of this reduced density matrix contains the subspace spanned by the 
solutions of 

d-l 

E =0. (11) 

ji,---,jk=0 

Let the null space of the reduced density matrix of k adjacent sites, denoted by 
Afe, be spanned by the orthogonal vectors Icq), (a = 1, • • • , s > d'^ — D"^). Then we 
can construct the local hamiltonian acting on k consecutive sites as 

s 

h:=Y^ Xa\ea){ea\, (12) 

a=l 
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where Aq,'s are positive constants. These constants together with the parameters of 
the vectors \ea) inherited from those of the original matrices Aj, determine the total 
number of coupling constants of the Hamiltonian. If we call the embedding of this 
local Hamiltonian into the sites / to / + /c by hi^i+k then the full Hamiltonian on the 
chain is written as 

N 

H = Y,hi,i+k. (13) 

1=1 

The state |\1') is then a ground state of this hamiltonian with vanishing energy. See 
[2T] for a more detailed discussion of the above points. 

A Hamiltonian derived as above does not have any particular symmetry. Indeed 
the above class include all types of Hamiltonians which have the matrix product 
state as their ground state. A subclass of these Hamiltonians however do have the 
symmetry of the ground state. Consider equation (llip . multiplying both sides of this 
equation by U~^{R) from left and U{R) from right, and using d?]), we find that if Ci...j 
is a solution of (llip . then Ri^i' ■ ■ ■ RjjiCi'^...j' is also a solution of the same equation, 
that is: 

Ri,i' ■ ■ ■ Rj,j'Ci' = ^{R)Ci,---j- (14) 

This means that the null space of the reduced density matrix is an invariant subspace 
under the action of the symmetry group R®^ . Thus the null vectors Icq,) transform 
into each other under the action of the reducible representation R®^. Such vectors can 
be classified into multiplets such that each multiplet transforms under one irreducible 
representation of the group R. Let the states transforming under the irreducible repre- 
sentation D^^ of the group, be denoted by |ep. Then the operators h'^ = \ 
is a scalar under the action of the group, that is 

[Z)'^,/i^] =0. (15) 

Hence to ensure the symmetry of the local Hamiltonian we write it as 

/i,,+fc = 5]A^/i^ = ^A^|e^)(e^|, (16) 

where the number of free couplings is equal to the number of multiplets which 
span the null space A^. 

3 A new representation for spherical tensors of arbitrary 
rank 

We are now equipped with generalities about matrix product states and their sym- 
metry properties. In this section we specialize the above discussion to construction 
of spin-s MPS invariant under rotation in spin space. For such a chain we take local 
Hilbert space to be spanned by the d = 2s-\- \ states of a spin-s particle, i.e. the states 
m),m = —s, ■ ■ ■ s}. Let us denote the D dimensional matrix assigned to the local 
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configuration m by Ag^m- Rotational symmetry in the spin space now demands that 
the matrices Ag^m form an irreducible tensor operator of rank s in the space of D 
dimensional square matrices. In view of Q, we should find 2s + 1 matrices Ag^j^ such 
that the following relations are satisfied 

[L+,Ag^rn] = \J s{s + 1) - m(m + l)^s,m+l 

[L^,As^rn] = \l s{s + 1) - m(m - (17) 

where L^, L-|_ and L_ are the D dimensional representations (not necessarily 
irreducible) of the Lie algebra of su{2): 

[L,,L±] = ±L± L_] = 2L,. (18) 

Remark: For simplicity, we will use the notation \m) and A^ instead of \s,m) 
and Ag^rn respectively, when the label s is clear from the context. 

It is crucial to note that it is not always possible to find tensor operators of a given 
rank for a given D dimensional representation. For example while there is tensor of 
rank one, in two dimensions, given by 

Ai = -V2a+, Ao = a„ A_i = V2a_, (19) 

leading to the spin-1 AKLT model [5], with parent Hamiltonian ([2]), there is no rank 
1/2 tensor operator in 2 dimensions. By this we mean that if we take D = 2, then 
there is no non-zero solution for the following system of matrix equations 

[L,,A±] = ±^A± 

[L+,A+] = 0, [L+,^_] = A+, 

[L_,A+] = A_, [L_,A_] = 0. (20) 

Therefore the first task for construction of rotationally invariant matrix product 
states for quantum spin chains or quantum ladders is to have a compact expression 
for spherical tensors of arbitrary rank. 

A possible procedure for obtaining spherical tensors of integer rank is to take two 
low-rank (possibly identical) tensors and decompose their ordinary or tensor product 
by the Clebsh-Gordon series to obtain irreducible tensors of higher rank. In fact if 
Ag^rn and ^s',m' are two spherical tensors, then one can form the product Ag^mAg'^m' 
(if their dimensions are the same) or Ag^m ^s',m' (otherwise) and decompose the 
products by using the Clebsh-Gordon coefficients to obtain spherical tensors of higher 
rank. For example take the AKLT tensor of rank one. Ordinary multiplication of this 
tensor, does not give a tensor of rank 2, since cr^ = 0, however its tensor multiplica- 
tion gives a tensor of rank two of dimension 4, i.e. ^2,2 = 2cr+ ® (T+, etc. In this way 
the product of two rank-1 tensors can be decomposed to give a rank-2, a rank-1 and 
rank-0 tensor. The obtained tensors can again be multiplied with other tensors and 
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decomposed to obtain tensors of even higher rank. This procedure however has several 
drawbacks: first the dimensions of the matrices wih grow very fast as we increase the 
rank of tensors, second it requires multiple use of Clebsh-Gordon coefficients which 
makes the final expression of the tensors, especially for high-rank tensors, quite cum- 
bersome and not useful. Another useful procedure, is to invoke the Wigner-Eckart 
theorem which decomposes the matrix elements of any spherical tensor in the angu- 
lar momentum basis, to an angular part, which is the Clebsh-Gordon coefficient and 
a reduced part, which essentially defines the tensor. However this procedure does 
not always lead to a compact notation for the tensor operators themselves and the 
multiplication of such tensors requires heavy use of Glebsh-Gordon coefficients. In 
this paper we introduce a compact and transparent formula for spherical tensors of 
rank s, for s integer or half integer, and use it to construct matrix product states for 
spin chains. For rank-s tensors the dimensions of the matrices are 2s + 2, thus the 
dimension grows linearly with rank. 

To construct the spherical rank-s tensor, let us take the orthonormal basis {|s, m), m = 
—s ■ ■ ■ s} of the spin s representation and augment it by the single state |0), of the 
spin representation 

(0|0) = 1, (0|s,m) = (s,m|0) = 0. (21) 

On this larger space, the following is the reducible s © representation of angular 
momentum algebra: 

s 

Lz = ^ m|s,m)(s,m| 

m=—s 

I/+ = ^ \J s{s + 1) — m{m + l)|s, m + 1) (s, m| 

m=—s 

L_ = ^2 ^-3(5 + 1) — m{m — l)\s , m — 1) {s , m\. (22) 

■m=— s 

Now it is readily verified that in this 2s -|- 2 dimensional space, the following 
matrices form an irreducible rank-s spherical tensor, that is they satisfy the relations 

As,m ■■= \s,m){ 0| + (-l)^-™|0)(s,-m| (23) 

where — s < m < s. 

It is important to note that the rank of these tensors can be integer or half integer. 
Such operators transform as an irreducible rank s tensor in the space which carries 
the reducible representation s © 0. 

Note: One can define the tensors more generally as 

As^m ■■= a\s,m){ 0| + 6(-l)'^|0)(s,-m| 

where a and b are arbitrary numbers, however these tensors are equivalent to the 
previous ones in the sense that they reduce to them by a suitable unitary transforma- 
tion. The factor (—1)'^ is inserted in the definition to ensure that no complex number 
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Figure 3: Color Online. The fully dimerized spin-s states in equation (I26p . 

enters the expression for half- integer ranks. 

While there are many representations for spherical tensors of different ranks, and 
these have been used in different works to construct various examples of invariant MPS 
dSl IZl dSl [m [211 ESI [26], to our knowledge the representation ([23]) is introduced 
for the first time. In the sequel we show that this representation is very general, 
in the sense that we can use it to find MP representations for all the variety of 
AKLT states, including the Majumdar-Ghosh or fully dimerized states of arbitrary 
spin, the partially dimerized states, and also the various states which are found from 
these partially dimerized states by different types of projection. Even more, one can 
construct other states not listed in the original AKLT papers, these are the symmetry 
breaking states. 

4 The spin-s fully dimerized or Majumdar-Ghosh states 

Using the definition of Am we find: 

N 

AmiAm2 ' ' ' Am2N — "'^'"^ '5m2i_i,-m2i |0) (0| 

4 = 1 

+ ..v (24) 

1=1 

Taking the trace we find 

N N 

Tr{Ami ■ ■ ■ Am2N) = "^^'~^^rn2i-i,-m2i + H^"''')^ '"^''^m2i,-m2i+i -(25) 

i=l i=l 

Inserting this into ([3][4]) the final simple form of the matrix product ground state is 
obtained as 

IV') = I0l) + 102) = \S)l2\S)34 ■ ■ ■ \S)2N~1,2N + \S)23\S)45 ' ' ' \S)2N,l, (26) 

where the singlet states \S) are given by 

\S) = ^=i== J2{-ir+^\m, -m). (27) 
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Figure 4: Color Online. 2s valence bonds (spin 1/2 singlets) in the AKLT con- 
struction, are replaced in our work, with a single spin-s valence bond, which is MPS 
representable by the matrix (j28p . 

Note that \S) is a singlet state, i.e. Lz\S) = L+\S) = L-\S) = 0. Thus {(pi) is a 
juxtaposition of spin-s dimers on sites (1, 2), (3, 4), • • • (2A^ — 1, 2A^) and \(j)2) is a one- 
site translation of |0i), i.e. a collection of spin-s dimers on sites (2, 3), (4, 5), • • • (2A^, 1), 
figure dS]). 

5 Other types of dimerized states 

A general dimerized state is one which is shown in figure (jH-a), where each line stands 
for a spin-1/2 dimer. The numbers of dimers are 2s and 2s' respectively. In our 
representation, we replace 2s spin-1/2 dimers with a single spin-s dimer, as in figure 
([D-b). Such a state has simple MPS representation, in the form 

\ip) = tr{Am^Am2Bjn[Bj^,^AmsAmiB^,^B^>^ ■ ■ ■)\mim2m[m2rn3rn4rn'^rn'4^ • • •) , 

(28) 

where the matrices {^4^} and {B^'} are embedding of the rank-s and rank-s/ tensors 
(1231) into a representation spanned by the vectors {|0), |s,m), |s',m')}, i.e. the direct 
sum representation s © s' © 0. In fact it is readily found that with 

As,m = Am = |s,m)(0| + (-l)^-'"|0)(s,-m|, 
^',m'^i?m' = |s',m')(0| + (-l)^'-"^'|0)(s',-m'|, (29) 

we have 

Am,A^,B^^'^B^'^ = (-l)'^-"^Mmi,-m2(-l)''""'^'^K.— (30) 
which readily yields the following partially dimerized form for the state ([28]): 

\i^) = \S)l2\S')34\S)56\S')78--- (31) 

where 15") and |«S") are respectively spin-s and spin-s' singlets defined in ()27p . 

One can construct even more general states, i.e. the symmetry breaking states of 
the form shown in figure ([5]) where the dimers are interspaced by spins which align 
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Figure 5: Color Online. The tensors (j23p can also be used to construct symmetry- 
breaking matrix product states. . 

in a particular direction. Consider the state 

= triA^^ A^^ C,„^ • • •) I mi , m2 , ms , m4 , ms , me , • • •> (32) 

where Am is of the form (j23p and Cm = am|0)(0|, in which a^'s {—s < m < 
s) are arbitrary complex numbers. Then the MPS represents a symmetry break- 
ing state shown in figure ([5|), where spins, 3,6, ••• are aligned in the state \a) := 
J2m=-s ^m\s,m) and the rest of the sites are dimerized. A suitable projection of 
these states, gives symmetry-breaking non-dimerized states [32] . 

6 The AKLT/VBS states 

In the AKLT models, one can use the fully dimerized states and project them to states 
which are called VBS states. While the parent Hamiltonian of the fully dimerized 
states has an interaction range of 3 sites, the VBS states which are obtained by 
projection have parent Hamiltonians with interaction range of 2 sites. The method 
is explained in figure ([6]), where we use a single spin-s dimer to replace 2s spin- 1/2 
dimers in the original method of AKLT. 

The lower state is obtained by projecting the states inside each bulb in the upper 
chain onto the symmetrized spin sector with total spin s + s' . It is now obvious how 
the parent Hamiltonian of the lower chain, the Hamiltonian which has this state as 
its ground state, should be constructed. Consider the first bond in figure ([6]) whose 
local Hamiltonian is denoted by hi. Due to the 2s' singlets, between the two bulbs, 
here we are only summing over two spin s states, instead of the apparent two spin 
s and two spin-s' states. Hence all the projectors Pj, with 2s + 1 < j < 2s + 2s', 
annihilate this bond, i.e. the local Hamiltonian hi, can be constructed as a linear 
superposition of all the above projectors with positive coefficients. By the same 
reasoning the local Hamiltonian /i2 can be a linear superposition of all projectors P^, 
with 2s' + 1 < A; < 2s -|- 2s'. Thus to construct a translation- invariant Hamiltonian, 
the parent Hamiltonian of the lower state can be constructed as 

H = Y,h,^,+i, (33) 

i 

where 

2s+2s' 

h= ^iPj^ (34) 

j=max{2s ,2s')+l 

where Pj's are the projectors on spin j sector of two sites and Xj are positive coeffi- 
cients. 
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Figure 6: Color Online. The projection method: From a partially or fully dimerized 
state, whose parent Hamiltonian has interaction range 3, an AKLT/VBS state is 
constructed whose parent Hamiltonian has interaction range 2. The form of the 
Hamiltonian is given in equations ()33p and ()34p . 

Note that the state on the lower chain is no longer dimerized, i.e. spins which 
are further apart than one lattice spacing, are correlated. Needless to say, the pro- 
jection method, although elegant in principle, is not suitable for calculation of many 
properties of the state. Having a matrix product representation for this state, turns 
all calculations into a straightforward and handy procedure. In this section we show 
that the irreducible tensors introduced in section ([3]), provides a MP representation 
for these states in a very simple way. 

The starting point of our procedure is however not to use 2s and 2s' spin-1/2 
dimers as in figure (jH-a), rather we use equivalently one spin-s and one spin-s' singlets 
as in figure (jH-b), for which we have already a MP representation. The spin-s and 
spin-s' dimers come from rank-s and rank-s' tensors (|23p . We multiply and symmetrize 
these two tensors to obtain a new tensor whose highest component is given by 

^s+s', s+s' — ^s,s^s',s' "I" ^s',s'^s,s- (35) 

From the explicit form of the tensors in ()23p . one sees that, 

Vs+s',s+s' ■= \s,s){s',-s'\ + |s',s')(s,-s|. (36) 

Note that this tensor lives in the 2s+2s'-|-2 dimensional space spanned by independent 
vectors {|s,m), \s',m'), —s <m<s, —s' < m' < s'}. 
It is readily verified that 

[L+,Vs+s',s+s'] = 0. (37) 

Therefore Vs^s',s+s' is indeed the highest-weight component of a spherical tensor of 
rank s -|- s'. Other components are obtained by successive commutations with L_. 
For example, we have 

Vs+s',s+s'-l ■■= J^—i{\s,S-l){s',-s'\-\s',s'){s,-S + l\) 
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The new spherical tensors have the interesting property that they lead to a non-empty 
null space A2. In fact it can be verified that these tensors have a peculiar fusion 
rule (decomposition of the product into irreducible representations), which exactly 
matches the fusion rule of the original 2s and 2s' valence bonds in a symmetric 
way. In the present formalism, this symmetry causes the final local Hamiltonian h to 
contain projectors common to both hi and /12 (figure ^) in the AKLT construction. 
Using the notation Vg to denote the whole multiplet Vs^rm ~s < m < s, the fusion 
rule of our tensors is 

m,ax{2s,2s') 

Vs+s'CSV,+,> = CjVj. (39) 

j=0 

Thus the multiplets Vj with max{2s,2s') + l < j < 2(s + s') are absent, i.e. identically 
vanish, in the decomposition of the left hand side tensors. In the language of matrix 
product formalism, section ([2]), this means that the null-space of the two-site density 
matrix, contains the multiplet of states which transform as spin j representations with 
max{2s,2s') + 1 < i < 2(s + s'). Therefore the local Hamiltonian annihilating the 
dimerized state, can be constructed from the projectors to these multiplets, namely 

h= ^^P^ (40) 

j=max{2s ,2s')+l 

where Pj's are projectors on spin j and Aj are positive coefficients. It requires tedious 
and lengthy calculations which may not be illuminating to prove (I39p in general. 
Instead we will give an idea of the proof by way of examples. First of all, it is readily 
seen from (f36]l that s, s' / 0, 

^s+s',s+s' — 0) 

but this is the top state of the multiplet V2s+2s' ^-nd hence this multiplet is absent in 
the right hand side of (139[) . In the same way one can also show from (]36p and (j38p 
that the top state of the multiplet V2s+2s'-i is zero. This pattern repeats until we 
arrive at the multiplet Knaa;(2s,2s')- We will give a more detailed and concrete example 
in section ([Tj). 

7 Examples 

Up until now we have been able to use our spherical tensors ()23p . in a uniform manner 
to construct all the variety of valence-bond states in the AKLT constrution. In this 
section, we will provide a few concrete examples. 

7.1 Properties of fully dimerized or spin-s Majumdar-Ghosh states 

First we calculate the normalization of fully dimerized states |^'^). The basic tool 
which we use is the following easily verified equation between the singlets, where 
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1, 2, 3 and 4 are any four different and not necessarily adjacent sites: 

^.^{S\S),,\S),, = ^-^\S),„ (41) 

where 

1 s = integer 

— 1 s = half integer 

This relation which we will use repeatedly in the following is depicted graphically 
in figure ([7]). Here a bulb around two sites 2,3 means that it has been multiplied 
from the left by a singlet 2,3 ('S'l- Repeatedly using equation (|1T|) or the graph ([7]), as 
in figure dSD, will give 



(42) 



(2s + l)^^i 
from which we obtain the normalization 

In order to find the correlations we use the following equations, 

^^{S\si-S2\S),, = -s{s + l), (44) 

and 

s(s + IV^ 

^^{S\s,.S2\S),,\SU = -s{s + l) ,AS\S),,\S),, = - ^J^^^ \S),„ (45) 

which readily gives 

((Ai|si -sal^i) = -s(s + l), ((/'2|si •S2|(A2) = 0. (46) 

Again the cross-product terms is calculated with the help of graph ([5]) , 

(01|S1-S2|02) = -g(g + l)(0l|'A2) = -Kg+l) (^^^-^^jV-l - (47) 

Putting these together we find the final form of the correlation functions: 
(^±|si • sal^^) s{s + 1) (2s + 1)^-1 ± 2e^ 



2 (2s + 1)^-1 ± ' 



(48) 



To construct the parent Hamiltonian of such states, we use (lllh and find the null- 
space of the reduced density matrices of three consecutive sites A3 , (two consecutive 
sites have no non-trivial null space in this model). From ()23p we have 

AmAm'A^,, = i-iy-"''6m',-m"H{0\ + 5„,^„„(-i)2«-™--" |o) . (49) 
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Figure 7: Color Online. The graphical representation of equation ()4ip . A bulb around 
2 and 3, means that the state has been multiplied from the left by 2,3(5*1, leaving the 
singlet on the right hand side. Note that the labels 1, 2, 3 and 4, denote any four 
sites. 




Figure 8: Color Online. The repeated use of the graph ([7]) for calculating the cross 
product terms {(j)i\(f)2)- 

To find the null space A3, we need to solve the matrix equation 

Cm,m' ,m" ■^m-^m.' -^ra" — 0) (50) 

m,m' ,m" 

which yields the following conditions: 

m 

5^(-l)™C™,_^,„, = Vm'. (51) 

m 

These conditions can be re-expressed in a more useful form, namely the null space 
A3 is spanned by vectors of the form 

\x) = XI Cm,m',m"\m,m' ,m") (52) 

m,m' ,m," 

which are perpendicular to the state (j27p i.e. 

i2{S\x) = 0, 2z{S\x) = 0, (53) 

where the subscripts indicate the embedding of {S\ into the local spaces of three con- 
secutive spins. Note that the factor \/2s + 1 has been inserted so that the state |S) be 
normalized. We will later use these equations to clarify the form of the Hamiltonian, 
but first let us derive an explicit form for the ground state. 

One is tempted to ask if or |02) are ground states separately. The answer is 
positive. To see this, note that the Hamiltonian is written in the form 

2N 

H = hk,k+i,k+2, (54) 
k=l 
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where h is the sum of projectors on the null space A3, i.e. 



a 



I Cq.) (Cq 



G Aa 



(55) 



Here {Icq,)} is a basis for A3 and from we know that (eQ,|5)i2 = (eQ.|5)23 V a. 
This implies that H\(pi) = H\(j)2) = 0. Each of the dimerized states and |(/'2), 
break the translational symmetry of the Hamiltonian. Finally let us also derive the 
parent Hamiltonians for the simplest cases, namely spin 1/2 which is well-known and 
spin-1 Majumdar-Ghosh states. 

a: The parent Hamiltonian for Spin- 1/2 dimerized state 

Using the standard notation |+) := |i, i), |— ) := |i, ^) we order the states of 
auxiliary space, as {|+), |0), | — )}. Then from (p3]) we have 



= |+)(0| + |0)(-| = 

A_ = |_)(o|-|0)(+| = 
which transforms as a rank 1/2 tensor with the generators given by 




(56) 
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(57) 



The singlet states are |5) = "^(1+) ~) ~ l~) +))• To find the parent Hamiltonian, we 
should solve equation ([5T]) . or what is the same thing, find states |x) G (81 (81 
such that i2('S'|x) =23 {S\x) = 0. It is readily found that there are four such states: 



|ei) = 1+,+,+), 

|e2) = -^(|+,+, -) + !+,-, +) + !-, +,+)), 

lea) = -^(1-, -,+) + !-,+, -) + !+,-,-)), 

|e4) = I-,-,-). (58) 

The vectors |ei) form the spin | multiplet, and if they come with the same coefficients 
in h in (j55p . the resulting Hamiltonian will be a scalar. It is known [221 [23] that in 
this case the parent Hamiltonian will be the Majumdar-Ghosh Hamiltonian, namely 
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b: The parent Hamiltonian for Spin-1 fully dimerized state 

Using the abbreviated notation |1) := |1,1), |0) := |1,0), |1) := |1,— 1), we have 
as the singlet state in \S) G , 

|S) = -^(|1,T)-|0,0) + |T,1)). (59) 

In order to find the null space A3, we note that due to su{2) symmetry, equation p^ . 
the basis vectors of A3 can be grouped into multiplets which transform irreducibly 
under su{2). These multiplets come from the decomposition of 1(8) I®! representation, 
which decomposes as 

1 (g) 1 (g) 1 = 3 e 2 © 2' e 1 © 1' © 1" © 0. (60) 

However not all the above multiplets belong to A3. In order to determine those which 
are, we should check the conditions ()5ip . It is sufficient to check these conditions only 
for the top state of each multiplet, since symmetry guarantees that the other states 
are present in A3, once the top state is present. With this insight we readily find the 
multiplets with the following top states are present in A3: 



ei) 


:= 1*3) = 


|1,1,1) 


62) 


:= N = 


|1,1,0) -|0,1,1) 


63) 


:= 1*2') = 


|1,1,0) -2|1,0,1) +|0,1,1) 


64) 


:= \ti) = 


|1,0,0) + |0,0,1) +3|1,1,1) + 3|1,1,1) -2|1,1,1) -4|0,1,0) 


e5> 


■■= \to) = 


|1, 1,0) - |1, 1,0) - |1,0, 1) + |1,0, 1) + |0, 1, 1) - |0, 1, 1), (61) 



where \tj) denotes the top state of the spin-j representation. One can verify that these 
are actually the top states by checking the equations Lz\tj) = j\tj), and L^\tj) = 
and also that they really belong to A3 by checking i2(5'|ej) =23 {S\ei) = 0. 

Having 5 different multiplets in the null space, means that the Hamiltonian has 
5 different couplings which can be tuned. Of course one of the couplings can be set 
to unity by a choice of energy scale. Let's call the projectors on the representation 
space j by Pj. Then the local Hamiltonian h will be 

h = XoPo + XiPi + A2P2 + X2P2' + A3P3. (62) 

Remark: It is important to note that the MPS formalism, gives a larger family 
of parent Hamiltonian than the original AKLT construction. In fact in the AKLT 
construction, the presence of projectors P3, P2 and P2' and Pq is automatic. However 
the presence of the new projector Pi is the result of the MPS formalism. 

The next step, which is not trivial, is to write the projectors in terms of local spin 
operators. The point is that on the decomposition (j60p only some of the representa- 
tions on the right hand side belong to A3. For those representations which occur with 
multiplicity one, we can easily find the expression of the corresponding projectors in 
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terms of local spin operators. Let us denote the sum of spin operators on three sites 
by S, i.e 

S := Sx + S2 + S3. 

The basis states of the representations on the right hand side of (j60p are such that 
they block-diagonalize the generators and hence the operator S • S. Let us denote 
the projectors on the totality of spin j representations by Vj, i.e. Vq = Pq, Vi := 
Pi + Py + Pyi , V2 = P2 + P2' ) and = P'i- Then we have the following system of 
equations (S • S)^ = T,f=oW + l))^^/> A; = 0, 1, 2, 3, or more explicitly, 

I = + + T'l + 7^0 
S-S = I2P3 + 6P2 + 2P1 
{S-Sf = 144^3 + 36^2 + 4Pi 

[s-Sf = 3437^3 + 1257^2 + 277^1. (63) 
Inverting the above equations we find 



72OP3 = 12S-S-8(S-S)2 + (S-S)3 

144^2 = -24S • S + 14(S • S)2 - (S • S)^ 

8OV1 = 72S • S - 18(S • S)2 + (S • S)3 

144^0 = 144 - 108S • S + 20(S • S)^ - (S • S)^. (64) 

A positive linear combination of the projectors 7^3, 7^2, and Vq gives a three parameter 
family of Hamiltonians. The projector Vi should be left out from this combination, 
since only one of the spin-1 representations belong to the null space A3. In general 
those representations which occur with multiplicity one, can always be expressed in 
terms of total spin operator S on three sites. However we can construct a more general 
family of Hamiltonians by calculating explicitly all the projectors in (j62p in terms of 
the most general set of independent three-body spin operators. A straightforward 
calculation gives the final form of the Hamiltonian (with the abbreviation 5*12 := 
81-82): 

2N 

H = ^^Jo + Jl Si^i+i + J2 + J3 Sf i^i + Ja Sli+2 + J5 5'j-|-i,i+2} 



i=l 



where 



J6{Si,i+2, {Si^i+l, 5j+i^i+2}} + J7 (5'i,i+l5'i_i+2'S't+l,t+2 + 5'j+l,i+25'j,j+25'i,i+l) , 



3 1 1 
Jo = -2Ao + -Ai + -(A2 + A20 + T-A3, 
5 3 15 

Jl = 2Ao--Ai--(A2 + A2') + T7A3, 
11 

J2 = -3Ao + — Ai + -(A2 + A2')- Y^A3, 

13 1 1 

J3 = 2Ao - — Ai - -(A2 - A2') + T7A3, 
20 6 15 
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Ji — Ao + — Ai + -(A2 - A2') + :j7;A3, 
20 ^ oO 

^5 = Ao — -(A2 + A2') + ttAs, 

1 11 
Je = -Ao +^Y0^i +^gA2' - 

J7 = Ao-^Ai-^(A2 + 5A20 + ^A3. (65) 

This Hamiltonian may seem complicated and not so interesting from the physical 
point of view. However we should note that it has effectively four adjustable parame- 
ters, (sine we can take Aq = 1) and by tuning these parameters this Hamiltonian may 
come close to physically simple and interesting models. For example if we take the 
parameters as follows: 

Ao:=l, Ai = 8 + ^A2, A'2 = 12 + 2A2, A3 = 18 + 4A2, (66) 

then the couplings J3, J4, Jg, and J7 all vanish and the Hamiltonian finds the following 
simple form, modulo additive and positive multiplicative constants 

2N 

H = '^Si ■ Sj+i + ASj • Sj+2 + K{Si ■ Sj+i , Si+i • 84+2}, (67) 

i=l 

where A := and K = ^x^f^- By taking A2 = 6, we can further set A = 0, 

and K = ^ and hence we can arrive at 

2N ^ 

-ff = ^ Sj • Sj+1 + ^{Sj • Sj+1 , Sj+1 • Si+2}, (68) 
i=i 

or by taking A2 very large, we can come arbitrarily close to the following Hamiltonian: 

2N ^ ^ 

-f^ = ^ Sj • Sj_|_i + -Si ■ Sj+2 + T7:{Sj • Sj+i , Si-i-i • 84+2}- (69) 



_^ 5 ^ 10 



7.2 Examples of AKLT/VBS states 

While the MPS representation may not be a necessity when dealing with fully dimer- 
ized states , such representation is invaluable when dealing with VBS states. 



Spin 3/2 VBS state 



As our last examples, we consider the MP representation of the spin 3/2 VBS state 
of the form shown in figure ([2]), which is obtained from a dimerized state with s = 1 
and s' = 1/2. From ()36p we see that the MP representation of such a chain is given by 
the following matrices, where we have abbreviated |l,m) — > |m), and |^,ib^) — > |ib) 

^3/.,3/. = |l)(-|+l+)(-l| 



19 



W = /f|o)(-| + v^|-)(-i|-/^|i)(+|-/fl+)(o| 

^3/.,-!/. = -vlio)(+i-/^i+)ai + /ii-i)(-i + /fi-)(oi 

^3/.,-3/. = -|-1>(+|+|-)(1|. (70) 

Note that we use equation ([36]) to find the highest- weight component of this tensor 
and the rest of the components are derived by action of L_. In a basis with the 
order |1, 1), |1, 0), |1, —1), 1+), |— ), the 5 dimensional vectors ^3/2,01 have the fohowing 
exphcit form: 



/ . . . . 1 \ 



3/2,3/2 



3/2,1/2 



1 

7! 



-V2 



■1 • \ 
. V2 



■ I 



(71) 



3/2,-1/2 



1 

7! 



-V2 



J 



3/2,-3/2 



V 1 



(72) 



From figure ([2]) and the discussion following it, we see that the translation- invariant 
parent Hamiltonian annihilating this state, should be constructed from the projector 
P3 onto spin-3 states. In the matrix product formalism, this means that the null-space 
of the two-site density matrix, should contain the multiplet of states which transform 
as spin 3 representation. In view of (|39p. this means that in the decomposition of 
quadratic product of tensors V3/2,m) into irreducible representations of su(2), the 
representation of spin-3 should not appear, i.e. the components of this tensor should 
identically vanish. This is indeed the case as one can see from (I72p that V^^^ ^ which 
is the component with highest weight of spin-3 representation vanishes. The other 
components vanish by symmetry. 

Equation (j39p generalizes this to arbitrary spins in a nice way which is exactly 
what we see in the valence bond picture of ([6]). Even more than that, it gives in one 
shot, the Hamiltonian which is common to both hi and /12 in figure 



The fusion rule of the tensors Vs 



As stated above the properties of valence bonds in the AKLT formalism are nicely 
captured in the fusion rule of the tensors Vg-, equation (|39p . Although the complete 
proof of ()39p is possible, we think it is not so illuminating. Instead we try to illustrate 
the idea by two simple example. Consider figure ([6]), with s = 3/2 and s' = 1. In our 
picture the VBS state obtained from projection is an MPS with auxiliary matrices 
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given by ^5/2 := V'3/2+1. Note that we use to denote the totality of matrices Vs,rm 
for all m. We want to show that 

y5/2'^V5/2 = Vo(BVi(BV2(BV3, (73) 

that is, we want to show that in the decomposition of the left hand side the tensors 
V4, and V5 do not appear, hence local Hamiltonian annihilating the state of the lower 
chain in figure ^ can be constructed from projectors P4, and P5. To prove this we 
need only show that the highest weight components of the tensors V5, and V4, in the 
decomposition of the left hand side of (j74p vanish identically. To this end, let us 
write the explicit form of the components of V5/2, obtained from the top component 
^5/2,5/2 by using ([36l) and applying the commutation relation Vs^m [L^,Vs^m+i]- 
Ignoring the numerical coefficients and signs in front of all states on both sides, which 
are irrelevant for the following proof, and using the shortened notation \s,m) — > \m) 
(i.e. |l/2) = |3/2, 1/2), | - 1) = |1, -1)), we have 

^5/2,5/2 = + 

%2,3/2 = l^)(-l| + l^)(0| + |0)(^| + |l)(^|. 

It is now easily seen that V^^2 5/2 ~ implying that the highest weight of the 
V5 vanish. Moreover we see that V5/2,5/2V5/2,3/2 = ^5/2,3/2^5/2,5/2 = 0, implying 
the highest weight of V4 also vanishes. This example corresponds to figure ([6]) with 
s = 3/2 and s' = 1 (or with 3 and 2 valence bonds in the AKLT construction). 
There is a very interesting point here which we should mention. The point is that 
a spin 5/2 VBS state can also be constructed in the same way as in figure ([6]) with 
s = 2 and s' = 1/2 or as in the original picture, from partially dimerized states 
with different numbers, namely with 4 and 1 valence bonds. Here we expect that the 
local Hamiltonian which is used in the construction of translation-invariant state be 
constructed only from the projector P5. This is nicely captured in the fusion rule of 
our tensors, W5/2 := ^^2+1/2, which is 

Wr,/2^W5/2 = W4eW3®W2®WieWo. (74) 

In fact we have (again ignoring numerical coefficients on both sides), and with the 
same type of shortened notation as in the previous example, 

^5/2,5/2 = |2)(-l/2|+|l/2)(-2|, 

^5/2,3/2 = |l)(-l/2| + |2,)(l/2| + |-l/2)(-2| + |l/2)(-l|. 

It is now seen that while the top state of is zero, the top state of that 
is ^5/2,5/2^5/2,3/2 - ^5/2,3/2^5/2,5/2 is non- vanishing, proving the fusion rule ([39]). 
This argument can be generalized to the arbitrary spins s and s' , although the proof 
will not be more illuminating than the example given above. 
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8 Conclusion 



The main emphasis of this paper has been on the rotational symmetry properties of 
matrix product states. To this end we have constructed a simple representation of 
spherical tensors of arbitrary integer or half integer rank. A spherical tensor of rank 
s is represented in a 2s + 2 dimensional space, hence the dimension of space, increases 
only linearly with the rank of the tensor. The introduction of these tensors have 
made possible a unified approach toward fully dimerized, and partially dimerized 
or AKLT/VBS states. In this way we have been able to find a matrix product 
representation for all the variety of valence bond states introduced in the original 
AKLT paper. Having such a matrix product representation makes the calculation 
of many properties of such states, specially the non-dimerized states quite easy and 
straightforward. Moreover a MPS representation is more powerful, since it will give a 
larger family of Hamiltonians compared with the AKLT construction, since it allows 
to include more projectors in the local Hamiltonian. This will then lead to more 
flexibility in approximating realistic interactions with parent Hamiltonians of matrix 
product states. We have demonstrated this for a spin-1 family of Hamiltonians with 
nearest and next-nearest neighbor interactions. Finally we should remind that the 
above constructions can be generalized to other symmetry groups like su{n). At least 
for a self-conjugate representation of su{n), whose weight diagram is symmetric under 
reflection, then we can define tensor operators in exactly the same way as in equation 
([23]) . namely: 

:= |m)(0| + |0)(-m|, (75) 

whee |m) is the n — 1 dimensional weight vector of that representation. Such a MPS 
representation may be useful for example in recent considerations of AKLT models as 
in [33\ [M] where su{2) valence bonds have been replaced with su(n) valence bonds, 
or in [M], where trimmer ground states with su{3) symmetry have been studied. 
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